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1 Introduction 

The accepted theoretical model of elementary particles and their interactions - the Stan- 
dard Model (SM) - has been very successful at explaining the experimental data available 
up until now. It describes all the observed particles - leptons and hadrons - that interact 
by electromagnetic, weak and strong forces, and the predictions of the model have been 
tested in the particle collider experiments up to the energy scale of around 100 GeV. How- 
ever, it is clear that some new effects should occur at higher energies. One of such effects 
is the highly anticipated Higgs scalar particle that is predicted by the SM and is yet to 
be found experimentally. Prom another point of view, it has been shown that the SM can 
not be valid for arbitrary high energies, therefore, new physics should be discovered at 
higher energies, possibly already at the scale of 1000 GeV. 

Theoretical physicists have been creating models that could predict this new physics 
at higher energies, one of the most realistic ones being the Minimal Supersymmetric 
Standard Model (MSSM). However, there is no experimental data available at higher 
energies to test the predictions of these new models (and to constrain the imagination 
of theoretical physicists), and that is probably the biggest problem today in the field of 
high-energy elementary particle physics. Of course, huge effort is being put into building 
particle colliders that would operate at higher energies: for example, the Large Hadron 
Collider (LHC) at CERN, which will reach energies of 14000 GeV, will start its operation 
within the coming year and the physics community is eagerly waiting for the results. It 
is clear, though, that the experimental data will be more and more difficult to attain at 
higher energies, and so we must look for alternatives. 

One of the alternatives comes from the very nature of quantum physics: due to the 
quantum fiuctuations of the energy of intermediate particles during a particle collision 
process, the physics at higher energies (e.g. 1000 GeV) affects the results of the exper- 
iments where the initial and final particles are at lower energies (e.g. 100 GeV). These 
effects, naturally, become weaker as the difference between the energy of the experiment 
and the energy of the new physics grows. But that brings up an important point: if we 
increase the precision of the experiments and the theoretical predictions, we can analyze 
these small effects, and thus analyze the physics at higher energies without actually hav- 
ing experiments at these energies! Such analyses are called precision measurements, and 
are becoming increasingly important, as the experiments at higher energies become more 
difficult to perform. 
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The theoretical predictions, that can be compared with experiments, are obtained by 
calculating amplitudes for the interaction processes. The amplitude can be represented 
by a sum of Feynman diagrams, while each Feynman diagram represents an algebraic 
expression, which is constructed using the Feynman rules of the theory. Now, the Feyn- 
man diagrams are classified according to the number of loops they contain into tree-level, 
one-loop level, two-loop level and so on, and each additional loop level gives smaller cor- 
rections to the amplitude than the previous one. A typical process in the SM would 
contain several diagrams at the tree-level, which can be quite easily calculated by hand, 
however already at one-loop level, which is essential for any precision calculation, the 
number of diagrams grows to several hundreds! Such calculation is practically impossi- 
ble to do by hand, therefore, it is necessary to develop programs that can generate and 
calculate Feynman diagrams automatically. 

There are already tools developed, that are capable of automatically generating and 
calculating Feynman diagrams at one-loop level in the SM and MSSM - the most notable 
being the FeynArts package [31 H] and GRACE [5] - but together with the development 
of such tools there is also a need for ways to check the validity of the computed results, 
in order to use them with any reliability. It turns out that one powerful tool to perform 
such checks arises from the procedure of gauge-fixing the theory, which is performed to 
deduce the Feynman rules from the Lagrangian of the theory. 

Gauge-fixing lets us freely choose the gauge in which to perform calculations, and 
this gauge choice can be parameterized by the so-called gauge-fixing parameters. The 
Feynman rules will then depend on the gauge choice, however, the final results have to 
be gauge- independent. What makes this interesting is that the gauge- invariance of the 
final result occurs through a combination of all the diagrams, while each one by itself is 
gauge-dependent, so if some diagram is missed, or some other mistake occurs, the final 
result will very likely become gauge-dependent. Therefore, if the result does appear to 
be gauge-invariant, it is a strong indication for the correctness of the result. The gauge- 
invariance check is also very suitable for automatic calculations, because it by itself can 
be performed automatically. If the Feynman rules with unspecified gauge parameters are 
used, the gauge parameters will carry over to the expression for the final result. The 
result can then be evaluated numerically with various choices for these parameters and it 
must not depend on their values. Naturally, the more gauge-fixing parameters we have 
available, the larger set of gauges we can explore, and the more stringent the test becomes. 
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The standard derivation of the vertices in the Standard Model uses linear gauge- 
fixing which results in gauges parameterized by 3 numbers {^w,^z,^a), that appear in 
the particle propagators. A natural extension of this procedure is to use a non-linear 
gauge-fixing which can result in a much wider set of available gauges with parameters 
that modify not the propagators but the vertices of the theory. The non-linear gauge- 
fixing of the Standard Model was presented in |8J, where it was used to simplify certain 
tree-level calculations. It was later recognized to be very useful for automatic calculations 
and it was fully implemented in GRACE [5], which was used then to perform one-loop 
level calculations in the Standard Model and test the results for the full non-linear gauge 
invariance. 

The purpose of this work is to implement the same class of non-linear gauge-fixing 
that was used in [HI [5] in the FeynArts/FormCalc package. The advantage of this package 
is that it is open-source and freely distributable, which means that the implementation, 
once tested to be correct, can be used by the whole physics community to check their 
calculations for gauge invariance. 

In the previous papers [6l [7] we have already explored the linear gauge-fixing in detail 
and tried to implement a subset of the non-linear gauges in FeynArts. In this final 
work we give a full presentation of the Standard Model, non-linear gauges, and their 
implementation in FeynArts. We start by discussing a Yang-Mills gauge theory with a 
broken symmetry, focusing on the gauge-fixing procedure. The results of this discussion 
are then applied to the specific case of the Standard Model with the non-linear gauge- 
fixing, which in the end yields a complete list of propagators and vertices of the theory. 
We then implement these new vertices in the FeynArts package and use it to carry out 
some automatic amplitude calculations and check the results for gauge invariance. 

2 Theoretical background 

In this section we present the general Yang-Mills (YM) Lagrangian with a broken sym- 
metry, which is the basis for most of the realistic quantum field theories, including the 
Standard Model, presented later in this work. The goal of this section is to familiarize 
ourselves with the YM Lagrangian and the gauge-fixing procedure, arriving at the final 
form of the Lagrangian, which can then be used in the next section for describing the 
Standard Model. 
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Note that the following discussion is fairly standard and closely follows the one found 
in [T], however, we present a more complete discussion of the linear gauge-fixing by con- 
sidering the most general linear case and consequently finding the largest possible set of 
linear gauge- fixing parameters. 

2.1 Yang-Mills Lagrangian and symmetry breaking 

We consider a gauge theory for a multiplet of real scalar fields 0i, transforming as some 
representation R of the gauge group G. The infinitesimal law of gauge transformation for 
the scalar fields 0j and the gauge fields is then [I]: 

5(t)i{x) = ia\x)t1.(t)j{x), (2.1) 
5A;{x) = -^d,a^{x) + r"'Al{x)a'{x), (2.2) 

where a"'{x) are the infinitesimal parameters of the transformation, tfj are the generators 
of the representation R, and /"^'^ are the structure constants of the gauge group G, 
meaning that the generators (in any representation) commute as: 

= z/^^'^f. (2.3) 

The parameter g is called coupling constant and can be chosen independently for each 
simple or U{1) subgroup of G. 

If the fields 0j are real, the generators t°- will have to be purely imaginary in order to 
make 6(j)i real. In this case it is convenient to define real generators 

T" = -it", (2.4) 

giving the transformation 

(50,(x) = -a'^(x)7;^.0,-(x), (2.5) 
and the commutation relationship 

pna^rj^b^^ jabcj.c_ ^2.6) 

In order to construct a Lagrangian, as usual in Yang-Mills theories, a covariant deriva- 
tive 

(D,0), = d,^, - igAltt^4>j = d.cp, + gAl^^cl)^. (2.7) 
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and a field strength tensor 



F;, = d,Al - d^A; + gr^^A^Al, (2.8) 



are defined. An invariant Lagrangian for this gauge theory is: 

4 ^ Ml"' ^ 2 



^ = -7(^;.)' + ^iD,4>,r - v{<p), (2.9) 



where V{(j)) is the symmetry-breaking potential. 

So far, in the Lagrangian above, the gauge fields are massless, but before proceeding 
to the quantization we have to include the effects of the spontaneous symmetry breaking, 
which occurs if the potential V{(j)) does not have a local minimum at = 0. A "normal" 
potential would have a minimum at = like V{(j)) = m^|0p, in which case the field 
fluctuates around 0, and the potential term gives the mass to the field. Now if we take 
V not to have a minimum at 0, like V{(f)) = a\(f)\'^ — the value of the field will "roll 
down" to the minimum of the potential at some (po, and fluctuations can only happen 
around that minimum value. In other words, the field will acquire a vacuum expectation 
value (VEV), 

{(l^oh = {<i>i). (2.10) 

which is saying that the field is non-zero even in the vacuum. 

In order for perturbation theory to work, we have to analyze fields that have zero 
vacuum expectation value, so we redefine (pi as 

(t)i{x) = (Poi + Xi{x), (2.11) 

and from now on treat Xi{^) as our field of interest. 



Plugging in the new definition (2.11) into expression (2.9) yields the Lagrangian in 



terms of Xi- Expanding it up to the quadratic terms we get: 

^2(A, x) = - \Al{-g^''d' + d^d'')Al + \{d,Xi)' 

+ g^^x^AlF\ + ]^g^F\F\AlA^^' - ^M,,xa„ (2-12) 
with constant matrices F": 

F\ = ^ct>,j, (2.13) 

and Mij the coefficient of the quadratic term in the Taylor series expansion of the potential 

r(0) at 0o: 
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Note that in the above Lagrangian the gauge bosons acquired a mass with the mass matrix 
(yf^F^jF^j, and the scalar fields have the mass matrix Mij. This is basically the famous 
Higgs mechanism, when the gauge bosons, initially massless, acquire a mass through a 
spontaneous symmetry breaking of the scalar field potential. 



In addition to the quadratic terms in (2.12), the full Lagrangian also includes various 
higher-order terms, which after quantization describe the interactions between the fields. 
For the remainder of this section we will not explicitly write those terms out, and instead 
focus on the quadratic terms and how they are affected by the gauge-fixing. The reason 
for that is that the quadratic terms define the propagators, and thus the particle content 
of the theory, and the higher-order terms is then treated as a perturbation of the quadratic 
Lagrangian. 



2.2 Gauge- fixing 

The Yang-Mills Lagrangian described above still can not be used for direct calculation 
of Feynman diagrams because it possesses a gauge symmetry, which makes the straight- 
forward path integrals badly divergent. The solution for this problem is the gauge- fixing 
procedure, which we describe here. As a result of it we will get an effective Lagrangian, 
which can already be used for calculating matrix elements. 



2.2.1 Faddeev- Popov procedure 

Quantization of a theory, using the path-integral method, proceeds by analyzing the 
quantity 

Z = j VAVxe'^^^'^'^'l (2.15) 

from which the propagators and interaction vertices can be read off. However, for a 
Lagrangian that is invariant under some gauge symmetry the straightforward analysis 



fails. The reason for this is basically because the integral in (2.15) runs over many field 
configurations that are gauge-equivalent, instead of counting each physical configuration 
just once. In order to avoid this problem, it is required to "fix the gauge", that is, to 
subject the fields to some constraint, which would prevent over-counting, and perform 
quantization under this constraint. 

Gauge-fixing is usually done by the Faddeev-Popov procedure, as a result of which Z 
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is expressed as [H |2] : 
Z = C ■ I VAVx exp 



det 



5a\y) 



a=0 



(2.16) 



Here G"[A, are arbitrary constraining function^ - there are as many of them as 
the number of gauge fields, and they are "constraining" in a sense that the condition 



= serves as a constraint on the fields in deriving (2.16). Then the new gauge- fixing 
Lagrangian term £gf is any function of the G", while C is an irrelevant normalization 
constant. Subscript a near a field means that the field is gauge-transformed with gauge 
transformation parameters a. 

The sum (£ + £gf) now looks as an effective Lagrangian, which we could just use in 
the calculation of diagrams, but there is still the strange determinant factor preventing 
us from that. Luckily, it can also be turned into an effective Lagrangian term as follows. 
It can be shown, that given any Hermitian matrix D"'''{x,y), its determinant can be 
evaluated by the field integral: 

det(Z,-(..,)) ^ jmVce^ (,J^.^y^i,)D'WM) . (2.17) 

where c"'{x) and c°'{x) are anticommuting fields. Using this relation, we can express (2.16) 
as follows: 

j d^x{C[A,x]+C^i{G-) + C^^[c,c,A,x]) , (2.18) 



Z = C ■ VAVxVcVc exp 



6a^ 



a=Q 



c\x). 



(2.19) 



Note that in the above equations compared to (2.17) one coordinate integral is missing, 
because the functional derivative is "local", in a sense that it contains a 5{x — y) 



factor, which cancels the integral. With Z expressed as (2.18) it is apparent that we can 
account for the determinant factor by introducing fictitious fields c and c, called ghosts, 



with their Lagrangian £gh as in (2.19) and then use the regular quantum field theory 
methods for calculations. Since the fields c where introduced as anticommuting numbers, 
the associated particles, when calculating Feynman, are to be treated in a way similar to 
fermions. 

In the end we see, that the consequence of the gauge-fixing procedure is the new 
effective Lagrangian 

>Ceff = C -\- Cg{ + Cgii (2.20) 



^Technically they are functional, because they can (and usually do) depend also on field derivatives 
at X 
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which has to be used instead of the original £ for the calculation of propagators and 
interaction vertices, also adding the new ghost particles to the theory. 

2.2.2 Linear gauge-fixing 



Let us now turn back to the Yang-Mills Lagrangian (2.12) and apply the gauge- fixing 
procedure, resulting in the effective Lagrangian. Here we will discuss linear gauge-fixing, 
meaning that the constraining functions are linear functions of the fields, and later we 
will consider the non-linear case. The usual choices for C and £gf for gauge theories is: 

G^{x) = l=[d^A^^{x) - ^9F%x^{x)], (2.21) 

£gf = (2.22) 

with an arbitrary gauge-fixing parameter C, - these are the so-called i?^ gauges. Such a 
choice is enough to do calculations, but it lets us explore only a one-dimensional space 
of gauges, while in the context of this work we are interested to have as much gauge 
choosing freedom as possible, that is, as many gauge-fixing parameters as possible. For 
that purpose here we will generalize the gauges by considering the most general linear 
constraint, constructed from gauge fields' four-divergences and scalar fields: 

= d,A-^ - K\xi, (2.23) 

with any real matrix K°'i. Note that any overall factors multiplying the constraints, or 
any linear redefinition of them (as G"" = L'^^G^) is irrelevant, because it can be absorbed 
in the definition of £gf . This is consistent with thinking of as constraints, because such 
redefinitions would actually yield the same constraints on the fields. The most general 
£gf we are then interested in is quadratic in G: 

£gf = (2.24) 

where L"-^ is any real square matrix and — ^ is for later convenience. First order of G 
should not appear in £gf, because the total derivative term d^A"'^ is irrelevant in the 
Lagrangian, and the first-order field terms Xi would shift the minimum of the potential, 
so we have to redefine the fields again. Then the higher order terms of G in Cgf would 
yield higher order field terms in Ces, which will be a matter of non-linear gauge-fixing. 
Another note: L"'' can be taken to be symmetrical, because the antisymmetric part of it 



cancels out in (2.24) and, therefore, is irrelevant. To summarize, our gauge-fixing choice 
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has now two matrices of parameters K and L (instead of a single number ^ ) . The goal now 
is to find all the relations these parameters have to satisfy, in order for the theory to be 
quantized properly, and to find the maximum number of really independent parameters. 
We now proceed with the calculation, using our definitions for and £gf. Plugging 



(2.23) into (2.24) we get: 



(2.25) 



which then together with (2.12) yields: 



1 



+ -Xii-d' - LK - M)i,xj 
- d^A'^^igF - LK)\xr 



(2.26) 



here we switched to matrix notation: = (F)''., K'^i = {Kf^, L"'' = {LY\ Mij = (M) 

We now enforce a condition on the Ces that the resulting propagators have to be 
diagonal in particle species. This is necessary in order for the fields in the Lagrangian 
to represent the physical particles (mass eigenstates) and for the perturbation theory to 
be valid. This diagonalization can be done by field redefinition, but it is much more 



convenient to use gauge-fixing to achieve that, when possible. The last term in (2.26) is 



a cross-term between A and x ^i-nd we have to get rid of it somehow. Obviously, we can 
use the gauge-fixing to do that easily by enforcing the first relation on our parameters: 



LK = gF, 



(2.27) 



or 



K = gL-'F, 



(2.28) 



leaving only L as the independent matrix. Here, taking an inverse, we assume that L is 
non-singular, so we will have to take a limit, if we want to consider some of the eigenvalues 



of L to be zero. Using the above expression for K, (2.26) becomes: 

Aff = - lA^i-g'^'-d^ + (1 - L)d''d^ - g^FF^g^T'K 



^2 g^F^L-^F-M),,x,- (2-29) 
One remaining ingredient for the effective Lagrangian is the ghost term £gh, which we 



calculate now. A variation in under a gauge transformation, using (2.23), (2.2) and 
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^ is: 

SG'' = d^iSA-f^) - K\{5x^) 
9 

= ^d^a" + r''^df,{A^''a') + K\F\a^ + K\T^^XjCi\ (2.30) 

Then the ghost Lagrangian is: 

£gh = c'^(x) - gK\F\ - gr' + A''^ix)d^) - gK\T^^Xj{x)) c\x), 

(2.31) 



where we multiphed (2.30) by an overall factor (— (?) - we can do that because it just 



changes the overall irrelevant factor C in (2.16). Let's again separate the quadratic term 
that gives the ghost propagators: 

£gh2 = c''{x){-d'^ - gKF^Y^c\x), (2.32) 

while the other two terms in £gh give ghost couplings to the gauge and scalar fields. 

To summarize, our final quadratic gauge-fixed Lagrangian, after enforcing the relation 
LK = gF and including ghosts is: 

CeS2 = - \A1{-g^''d^ + (1 - L)d^d^ - g^FF^gn^'Al 

+ \xi{-d^ - g'F^L-'F - M),,Xj (2.33) 

We can easily read off the propagators for the gauge, scalar and ghosts fields from this 
Lagrangian, but they will not necessarily be diagonal in the species space, and, therefore, 
the fields in this Lagrangian will not represent the physical particles. In order to proceed 
with finding the physical fields and their (diagonal) propagators, we have to know the 
F matrix, appearing in the Lagrangian, and so we will continue this analysis when we 
consider the Standard Model in particular. 



2.2.3 Non-linear gauge- fixing 

Linear constraints are the simplest way to perform gauge-fixing. Linear gauge-fixing is 
also essential because linear yields quadratic Cg{ and so it participates in defining 
physical fields and propagators. However, according to the Faddeev-Popov procedure 
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of gauge-fixing, can be any functions of the fields, and so in particular they can 
contain non-linear combinations of fields. Of course, we don't want to introduce non- 
renormalizability into the theory by the gauge-fixing procedure, so we have to keep the 
total mass dimension of the fields in £gf terms up to four. Additional requirement is 
that we would like to keep the quadratic part of £gf unchanged, since we have already 
determined how it should look. 

Assuming we keep £gf the same function of G", the modification that respects the 
requirements above is adding a quadratic term to G°^: 

= d^A^^' - K\xi + Gl,, (2.34) 

where G^j is the non-linear, that is, quadratic part. In this case the quadratic part of £gf 
will remain the same as before, but it will have terms of third and fourth order, that will 
give contributions to the interaction vertices. Terms of third and higher order in fields 
are not allowed in G", because vector and scalar fields have mass dimension 1, so that 
would lead to higher than four mass dimension in the Lagrangian. 

Apparently there is a lot of freedom for choosing the non-linear gauge-fixing terms 
G^i, which can lead to various modifications to the interaction vertices of the theory. We 
will explore those choices in detail in the following sections of this work, when we discuss 
the Standard Model specifically. 

To summarize, the full effective Lagrangian of the Yang-Mills theory is 

1 



= - -Al{-g>'^d^ + (1 - ^d^^d" - g^FF'^g^'^Y^A^ 
1 



, ^X^{-^'-9'F^L"^F-M),^XJ (2-35) 
+ e{-d^ - g^L-^FF^Y^c''+ 

with the quadratic terms written out explicitly. Note that the constant g and constant 
matrices F and M are the parameters of the specific theory, while the matrix L is an 
arbitrary matrix of gauge-fixing parameters - we will specify those constants in the next 
section for the case of the Standard Model. The terms of higher than quadratic order 
that describe the interactions are unspecified, but grouped in three categories: Ci is 



the interaction part of the original Lagrangian (2.9) containing terms of higher than 



quadratic order, are the interaction terms of the gauge-fixing Lagrangian arising from 
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the nonhnear constraints GjJi, and £ghi are the interaction terms involving ghosts which we 



have seen in (2.31), with some additional terms from the non-linear part of gauge-fixing. 

We will use this effective Lagrangian in the next section, applying the framework of 
this section to the specific case of the Standard Model. 



3 The Standard Model 

In the previous section we have described the Yang-Mills gauge theory with a broken 
symmetry in general, and showed how the procedure of gauge-fixing is performed, leading 
to a new effective Lagrangian. Here we will describe the Standard Model, also called the 
Glashow-Weinberg-Salam model, which is the accepted theory of electroweak interactions, 
and which is, in fact, a particular case of a Yang-Mills theory with a broken symmetry. 

The goal of this section is to give a theoretical background for amplitude calculations 
in the Standard Model and arrive at the full explicit expression for the Lagrangian. The 
most important section in the context of this work is the discussion of non-linear gauge- 
fixing, which gives us the new interaction terms, that we later implement in the FeynArts 
package. 



3.1 The boson content and the propagators 

Here we will work out the details of the quadratic part of the full effective Lagrangian 



(2.35), which will tell us the physical particles of the theory and their propagators. First 
we will determine the matrices F and M and then diagonalize the propagators, in order 
to find the physical fields. 

The electroweak theory goes along the general framework of the previous section, with 
the gauge group 

G = SU{2)xU{l), (3.1) 
with structure constants (assuming is the generator of U{1)): 

, e"^", a,b,ce |1,2,3}, 




The real scalar fields (pi are taken to transform as the real components of the usual two- 
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dimensional complex representation of SU (2) with generators: 



-, a £{1,2,3} 



1 



(3.3) 
(3.4) 



where cr" are the Pauli matrices and is the generator of U{1). We define the real 

components of the complex two-dimensional vector as: 



V2 I 04 + icf)^ 

then the real representation matrices T" = —it"- in 0' space look as: 



(3.5) 
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and their commutation relations are, naturally: 



(3.6) 



^rpa^rpb^^^abcrpc^ O, 6, C G { 1 , 2, 3} , 



(3.7) 
(3.8) 



Now, following the previous section, we consider the symmetry breaking of this group, 
which happens because the scalar fields are subjected to some potential V{(f)). A renor- 
malizable potential that obeys the gauge symmetry and doesn't have a local minimum at 
= 0, as required for symmetry breaking, is: 

■2 A 



n0) = -y(W + 4(W, 



(3.9) 



with some unknown parameters /i and A. As a result of it, acquires a vacuum expectation 
value (VEV) at the minimum of V: 



|0o| = V 



(3.10) 
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Note that all directions of are equivalent under the gauge transformation, so any choice 
of direction for 0o gives the same results, so we can arbitrarily choose it as: 



(3.11) 






\v J 

It is important to note that there is one unbroken symmetry direction: the combination 
of generators 

Q = T^ + T^ (3.12) 

leaves (po invariant, and, therefore, this transformation is still a symmetry of the theory, 
even after (p acquires a VEV. 



Now we calculate gF^'i according to (2.13) to be: 

/ ^ \ 
0^00 



\ 



(3.13) 



0^0 
-g' J 

here down the matrix is the a-dimension, and to the right is the i-dimension. Note that 
we introduced two different coupling constants, g and g', because, as mentioned earlier, 
they can be chosen independently for each semi-simple or U{1) subgroup. 



3.1.1 Gauge bosons 

Let us now continue with the task of diagonalizing propagators and first consider the 
gauge bosons, with their Lagrangian 

Ca = -^A^^i-g^'d^ + (1 - L)d^'d'' - g^FF^g^'^^Al (3.14) 

Because of the different spacetime tensor coefficients of L and of the mass term FF'^ each 
one has to be diagonalized separately. We first proceed with the mass term, because it is 
already given: 

/ ^2 \ 
^2 



g^FF"" 



V 

T 



g^ ~gg' 
\ -gg' g" J 



(3.15) 
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(3.16) 



The diagonalization of this term is done by field redefinition, so we want to express the 
mass term as: 

with unitary transformation U, and diagonal mass matrix Ma- We take: 

f \ 
^2 

g^ + g'^ 
\ 



Ma 



V 



(3.17) 



^71 

1 



U 



\ 









I 





i 








9 







g' 










g 



(3.18) 



which yields physical fields A' — UA and their masses as: 



Wr. 



V 



/3 



7y — A 



4 = 4"' 



(^4 - g'Al) mz = 



7'2_ 



V 5 + 



(3.19) 



niA — 0. 



This is the main result of the electroweak theory (independently of the gauge-fixing de- 
tails), that there are three massive gauge bosons, mediating the weak interaction, and 
one massless boson, which represents the unbroken symmetry and is nothing else but the 
photon. 

Before going on further let's define some new constants for later convenience. Consider 
replacing the two unknown coupling constants g and g' with other two constants: electric 
charge e and Weinberg angle 9yj defined as 



99' 



(3.20) 



with the inverse transformation: 



sin 0,, 



9 



cos 9,, 



(3.21) 
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It is not apparent yet, but it will be shown briefly that e as defined here is indeed the 
electromagnetic coupling constant. The definition of 9yj allows to rewrite Z and A more 
clearly as: 

Z^ = cosd^Al-smd^Al, (3.22) 
A^ = sin e^Al + cos e^A^, (3.23) 



and the masses as: 



mw = v , mz = v- — ^— . (3.24) 

2 sm 2 cos Uo„ sm &„, 



Also note an important relationship (or, alternatively, a definition of 6^): 

my/ = t^z cosOyj. (3.25) 



One more expression that will be very useful later is the covariant derivative from (2.7) 



D, = d,- ig{Ali' + 4? + Alt') - ig'Ali' (3.26) 

expressed in terms of the physical gauge fields. Note that the f here are generators and 
depend on the representation. Substituting A^ with the physical fields in the end we get: 

D^ = d^~ t^{W;:i^ + W;i-) - i-^Z^{P - sln^ e^q) - teA^q, (3.27) 

V 2 cos Ow 

with the new generators: 

i^={i^±iP), (3.28) 
q={P + i^). (3.29) 

We can see now that we have recovered the regular coupling of the field A^ with a coupling 
constant e and a charge given by the value of the generator q. 

Let's now finish analyzing gauge bosons. Some additional explanation is needed, why 
the are defined exactly like they are. Because there are two identical eigenvalues 
of the mass matrix, any unitary transformation of A^ and A'^ would lead to two energy 
eigenstates with the same mass. Therefore, additionally we want to make them eigen- 
states of the remaining symmetry generator Q, which with the identification of as the 
electromagnetic field, happens to be the electric charge operator. Since gauge fields form 
an adjoint representation of the gauge group G, we can calculate their transformation 
properties using commutation relations of the group operators. Specifically, if we define 
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the field operator, corresponding to field A" to be a linear superposition of generators 
(here we will use hats for operators to avoid confusion) : 

A = f^A"", (3.30) 

then the transformation of the field operator under a gauge transformation is given by 
the commutator 

M = -a"[f",i]. (3.31) 

We can use this relation to find the transformation of the fields W^. First we need the 
corresponding generators 2V± ■ Since the full field operator 

A = f''A'' = f^W^ + f-W-+f^Z + f'^A (3.32) 

is constructed from multiplying the fields with generators, the generators undergo an 
inverse transformation to that of the fields, therefore: 

(3.33) 

Now we can see that the commutators corresponding to the transformation of under 
the remaining symmetry: 



= —={f^±if'^). 
v2 



5f+ = -a[Q, T+l = -^[T^ + f\ + if 2] = --^(f 2 - if 1) = +iaf+, 

Sf- = -a[Q,f-] = -^[f^ + f\f^ - if] = -^{f^ + if^) = -iaf-. (3.34) 

indeed show that are eigenstates, or, to put it differently, transform irreducibly under 
the charge transformation. 

We can now put final constraints on the gauge-fixing parameter matrix L. With the 
gauge fields transformed to the physical basis as above, the Lagrangian for them becomes: 

jCa = - liA'y^i-g^^d^ + (1 - ULU^)d''d'' - MAg^"')^^A'^. (3.35) 

Since we can not use field redefinition anymore, we have to insist for the gauge-fixing term 

U = ULU^ (3.36) 

to be diagonal in this basis as well. This leaves us with four real parameters, which we 
write, in order to conform with the usual notation, as: 

/ i/e+ \ 
1/e- 
l/Cz 

\ 1/u J 



(3.37) 
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This gives us the original matrix: 



L = U^L'U 



—I 



2?+?- 







V 






















gg'(fe-€A) 



Since L can neither be complex nor asymmetrical, we conclude that 



giving the L matrices: 



\ 






V 




























\ 




^lu I 




99\(,z-U) 








and the other gauge-fixing parameter K: 



( 



K\ = gL-'F 



V 

2 



v 



















g^z 
-g'^z 






0/ 



(3.38) 



(3.39) 



(3.40) 



(3.41) 



(3.42) 



This is our final result for the parametrization of linear gauge-fixing. It leaves us with 
3 real parameters {C,W:^z,(,a), a slight generalization of usual gauges, which would 
correspond to taking all 3 parameters to be equal ^ = = ^z = ^a- 



Finally, the Lagrangian for gauge bosons in the physical basis (3.35) becomes explicitly 
diagonal: 

Ca = Y1 -lx;{-g^''d' + (1 - ^x)d^d'' - m\gnX., (3.43) 



X 



with X G {W^^Z^^A} and masses as in (3.19). The momentum-space propagators can 
now be easily computed by taking the inverse of the coefficient for the quadratic term in 
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the Lagrangian, and multiplying it by i. Switching to momentum space, we replace all 9^ 
by —ikfj^, which yields for gauge bosons: 

-i ( „., k^h}" 



/c^ — 



3.1.2 Goldstone bosons and Higgs 



Xe{W^,Z\A}. (3.44) 



Let's move on now to the quadratic scalar part of the Lagrangian (2.35): 
An important separation occurs in this Lagrangian as follows. It is a general theorem 



(3.45) 



(called Goldstone theorem) that the mass matrix of Xi defined in (2.14) satisfies 

(TyohMi, = F\Mi^ = 0, (3.46) 

for all a. It follows from the fact that V^(0) is gauge-invariant and it basically says 
that V{(j)) doesn't form a quadratic potential in the directions that 0o transforms to. 
Consequently, the only non-zero term allowed in M is: 



M44 = rn^j. 



We can confirm the above result using the explicit expression for V{(j)) from (3.9): 

dV{(t)) 



d(f)id(f)j 



00 



With 



vSjA this becomes: 



(-/i%- + A(0V)5.,+2A0V)| 



(w^A - /i^)5ij + 2t;^A5i4(5j4, 



and after substituting v = \/ /x^/ A: 

confirming the argument above with identification: 

TJlH = V2fi 



(3.47) 



(3.48) 



(3.49) 



(3.50) 



(3.51) 



On the other hand, the {g^F L ^F)ij term in (3.45) will be zero where i = 4 or j = 4, 



because of the zero column in F. Therefore, we can write as: 



1 



d'' - g^F^L-'F)mnXn, 



(3.52) 
(3.53) 
(3.54) 
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where m,n & {1,2,3} and we have renamed X4 = h. The Lagrangian separates and we 
get three massless fields Xm called Goldstone bosons (the mass-like term is proportional 
to L, which indicates that they have no physical mass), and one massive field h with 
an unknown mass rriH, which is called the Higgs boson. We see that the massive field is 
unaffected by the gauge-fixing procedure, and it already has a diagonal propagator 

i 



{hh) 



m 



(3.55) 



H 



It's worth noting here that the Higgs boson is a physical particle, which is an essential 
prediction of the Standard Model and which yet has to be discovered experimentally in 
particle colliders. On the other hand the three Goldstone bosons are not actually physical, 
that is, they can not appear in the initial or final state in amplitude calculations. Instead, 
they represent the missing degrees of freedom needed to make the three gauge bosons 
W^, Z massive from the initially massless vector bosons. This is in general the way that 
the Goldstone bosons function in theories with a broken gauge symmetry, and there will 
always be as many of them, as the number of massive gauge bosons. 

Of course, we still have to include the Gauge bosons as the intermediate particles in 
amplitude calculations, so we continue now with diagonalizing their remaining Lagrangian: 

1 



G 



with the mass matrix: 



/ .2 



9 iw 





w 



iXrn 








(3.56) 



(3.57) 



{g' + g")iz j 



We see that it is already diagonal in the current basis and the masses are ^/^mx, 
with first two Goldstone bosons corresponding to W and the third to . However, like 
with gauge bosons, since two mass eigenvalues are equal, the states should be chosen as 



eigenstates of charge operator, which, using (3.6), in the representation of Xm is: 

/ -1 \ 

Q = T^^T^= +1 . 

\ / 

We use the same redefinition for Goldstone boson fields as for gauge bosons: 

1 



(3.58) 



X 



± 



V2 



(Xi T 



(3.59) 
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which means that the states, similarly to generators of gauge bosons, undergo an inverse 
transformation : 

|x^> = ^(IXi)±^IX2)). (3.60) 
We can check now that these are indeed the charge eigenstates: 



|x^> = -aQ \X^) = -^i\X2) T I \xi)) = ±^^(IXi) ± ^ 1X2)) = ±m |x^> , (3.61) 



so the Lagrangian in the diagonal form looks as: 



(3.62) 



X 



with X e {H^ , Z }. That gives the following propagators for Goldstone bosons: 



{XxXx) 



A;2 - ^xm\ 



2 ' 



X e {W^,z"} 



(3.63) 



3.1.3 Ghosts 



Finally, let's look at the quadratic ghost Lagrangian piece in (2.35) 
It has the mass matrix: 



(3.64) 



V 



( 9\w 

g\- 



w 








\ 





(3.65) 



g^^z -gg'^z 

\ -gg'^z g'^^z J 

which, like the matrix for gauge bosons, is not diagonal. Notice, however, that it looks 
very similar to the gauge boson mass matrix, except for the ^ factors. Indeed, if we use 



the same unitary transformation (3.18) for the ghost fields, as we did for gauge bosons: 



(f/t) 



ba 



(3.66) 



then the transformed mass matrix becomes diagonal: 



^ g'^w 






V 



o\ 

g'^w 

ig' + g")^z 

y 



(3.67) 



24 



The masses for the ghost fields are then, similar as for Goldstone bosons, y^mx-, where 
mx is the mass of the corresponding gauge boson. However, unlike the case with the 
Goldstone bosons, there is also a massless ghost corresponding to the massless vector 
field A^j_. And so the diagonal Lagrangian for ghosts is then: 

= - ixm\)cx, (3.68) 



X 



with X G {W^ , Z ,A}, and the propagators are 



{cxcx) = _ ^^^y XG{l^^Z^A}. (3.69) 

3.2 Fermions 

Up to now we have only included scalar and gauge vector fields in our theory, because 
they are the crucial ingredients of a gauge theory with Higgs mechanism. In addition to 
that the Standard Model, of course, contains fermions which are the matter constituents, 
and so here we will discuss their Lagrangian. 

Let's start with a free massless two-component Weyl spinor tp, which can be either 
left-handed or right-handed. It has the Lagrangian: 

C = itp^tp. (3.70) 

In a gauge theory that extends naturally to 

C = i^p^ (3.71) 



with a covariant derivative as in (2.7): 



D^ = d^- igAlt\ (3.72) 

with generators dependent on the representation that a particular fermion belongs to. 
The only question then is how to group the fermions into representations of the SU{2) 
and what charges of the U{1) subgroup to assign to those groups. Note that the "charge" 
of U{1) is just the value of the generator t^, also conventionally called Y. 

We have the following observed fermions to fit into the Standard Model: electron e~, 
neutrino u^, up-quark u, down-quark d, and two copies of these particles with exactly 
the same properties but different masses called the second generation (/i~, z/^, c, s) and 
the third generation (r~, z/,-, t, 6). Note that all the particles here except for neutrinos 
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have both left-handed and right-handed components, but it appears that these compo- 
nents have to be treated as basically two distinct particles, for example and e^, only 
connected by a mass term, as will be shown briefly. Neutrinos have only the left-handed 
component, and we will skip the subscript L for them. 

The particles enumerated above arrange into representations as follows: the left- 
handed components form doublets of SU (2) 

called the lepton doublet and the quark doublet with U{1) charges 

Ye = -1/2, Yq = 1/6 (3.74) 

respectively. The right-handed components are singlets of SU(2): 

Gr, ur, dn, (3.75) 

with charges 

Fe = -1, K = 2/3, Yd = -1/3. (3.76) 

Remember that the components of the doublets have t^ = ±1/2, while for sing lets t^ = 0, 
and that the electric charge q = t"^ + Y, which gives the right charges = 0, qe = — 1, 
qu = 2/3, qd = —1/3. The other two generations look exactly the same, so we will not 
repeat them here, just keep in mind that all the calculations here apply equally to all 
three generations. 

We can now write the massless Lagrangian for all fermions in the first generation: 

C = iElPEl + lenlpeR + iQlPQl + iurPur + idapdR. (3.77) 



Using the expression for given in (3.27) and the values for the generators (t* = a"^ /2 
for doublets) we get explicitly: 

C =iEL^EL + ieR(^eR + iQL<^QL + iuR<^UR + idR^dR (3.78) 

=-L(z7,7^eL + ULl^dL), (3.79) 

J(^- =-L(ei7^zy, + dL^UL), (3.80) 

J^ = - e7^e + |m7'^m - ^d-f'^d, (3.81) 

[Ive^i^e - leil^eL + \ul1^ul - \dLl^dL - sin^ J^] . (3.82) 

COS f/^ii 
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The different fermion currents Ji^, J|^, J^, that were defined here each couple to a 
different gauge boson. Also note that the fields ^ m, c? without a subscript L or R denote 
the sum of both left (e.g. ei) and right (e^) components, and occur if the coupling is 
non-chiral. We can see that the electromagnetic current J^, as expected, is completely 
non-chiral, while the has a non-chiral part proportional to and a pure left-handed 
part. 

We have now the correct coupling of fermions to gauge bosons, but the fermions are 
still massless, which in reality is not true. The mass term for a fermion (e.g. electron) 
looks like 

£ = -rrieee = -meieieR + e/jez,), (3.83) 

but such term is not allowed in the gauge-invariant Lagrangian, because the left and 
right components belong to different representations of SU{2), thus their product is not 
invariant. In order then for fermions to acquire mass, we again use the Higgs mechanism 
- that is, we couple the fermions to the scalar fields (p and, after the symmetry breaking, 
the non-zero expectation value will provide the effective mass term. The allowed gauge- 
invariant coupling Lagrangian is [1]: 

£ = -XeELCpCR - XdQL<pdR - Xue'^'QLa^UR + h.c, (3.84) 

with unknown parameters Ag, A^, A^. The scalar field (f) here is the same two-component 



complex doublet under SU{2) as in (3.5): 



V2 \ (f)'^ + i(f)^ I V2 \v + h + ix3 



(3.85) 



In general, since such terms are not forbidden by any symmetry, there is no reason to ex- 
pect that they shouldn't be included, so their addition to the Lagrangian is fully justified. 
After expanding the Lagrangian we get: 

C = ^ (AeCe + Xddd + XuUu) [v + h) (3.86) 

V2 

+ i {Xei^LeR + XdUidR - XuURdL) 
+ i {-XeCRUL + XudiUR - XadRUi) 

- {XeeieR - XeCRCL + XddidR - X^dRdi - X^ulUr + XuUrUl) Xs, 



^not to be confused with the couphng constant e! 
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and the terms multiplied by v on the first line give masses: 

XuV XdV 

We can now freely take the masses mg, rriu, rrid as the unknown parameters instead of the 
couplings, and express the couphngs as 

Xi = = — . 3.88 

V y/2mw 



Combining (3.78) with (3.86) gives the full interacting Lagrangian for the fermions of the 



first generation: 

£/i =e{i(^ — me)e + vjyifp^^e + — ra^u + d{i(^ — rridjd (3.89) 

with fermion currents that couple to the Goldstone bosons and Higgs: 

J,+ =^jyLeR + ^UidR - ^URdi, (3.90) 

J- =- ^crVl + ^dLUR - ^dRUL, (3.91) 

j3 =JZ^e7^e + ^d-i^d - ^ufu, (3.92) 

Jh =^ee + ^dd + ^uu. (3.93) 

For the full fermion Lagrangian we then just need to sum over all three generations: 

3 

^/ = E^/- (3.94) 

i=l 

The only difference between the different generations are masses, which arises from the 
fact that the fermion-scalar coupling constants A can be different for each generation. 

There is actually an additional complication related to the inclusion of other gen- 
erations: the coupling to the scalars can actually mix generations, that is, there is a 
term X^lQLi'pdRj, where z,j are generation indices and is now a matrix. That effec- 
tively yields flavor mixing of quarks in the charge-changing weak interactions (through the 
so-called Cabibbo-Kobayashi-Maskawa mixing matrix), and allows a weak CP violation. 
Those effects are not, however, of our interest here, and we will assume A's are diagonal, 
which is a pretty good approximation. 
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Another note: the quarks, of course, come in SU{3) triplets, which is an unbroken 
gauge symmetry of the Standard Model, and so everywhere in the Lagrangian where there 
is a quark, it's actually a sum of three different color quarks. They have to be summed- 
over if a quark appears as an intermediate particle in amplitude calculations. There are 
then also gauge-bosons associated with the SU{3) called gluons and quark coupling, but 
we will not discuss that Lagrangian, as we will not be doing QCD calculations. 



3.3 Interactions 

Now we move on to the task of writing out all the remaining interactions in the Standard 



Model. Let's look back at the full effective Lagrangian for a Yang-Mills theory (2.35) and 
work on the three interaction terms Ci, Cgfi, £ghi in the following sections. The goal here 
is to have in the end a full explicit Lagrangian in terms of the physical fields, from which 
the interaction vertices can be read off directly. 



3.3.1 Ci - original interactions 

The term Ci contains vector-vector interactions arising from the dynamic terms of the 
gauge fields, vector-scalar interactions from the gauge coupling through the covariant 
derivative, and, finally, scalar-scalar interactions from the Higgs potential. To find the 



expressions for those terms we have to go back to the very original Lagrangian (2.9): 

^ = -l{F^,.r + l{D,<j^^y-V{<f>). (3.95) 

Writing it out more explicitly after the symmetry breaking it becomes: 

1 



C = --id,A:-d.A^^ + gnA''^A:r 

+ lid.x^ + 9AIF\ + gA^x,? (3.96) 



+ y(0oi + Xi)'- ^(</>oi + Xi)^- 

We can now collect the higher-than-quadratic terms, grouping them by the particle species 
that are involved in the interaction (e.g. C^w for 3-vector interactions and C^ss for vector- 
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vector-scalar-scalar interactions) : 



-^i ^VVVV ~l~ ^VVV ~l~ ^VVSS ~l~ ^VVS ~l~ ^VSS ~l~ ^SSSS ~l~ ^SSS") 


(3.97) 


2 


(3.98) 




(3.99) 




(3.100) 


'-'vvs — y ^ i-'-ijX.j-i 


(3.101) 




(3.102) 


^ssss ^ (XiXi) ) 


(3.103) 


= -\{(t>QiXi){XiXi)- 


(3.104) 



All the symbols used here (structure constants /, generators of the scalar fields T, vacuum 
expectation value matrix F — T^o) are defined in the earlier sections, so the only work 
that needs to be done here is to substitute all the definitions and the physical fields, and 
work out the final expressions. The computations are quite cumbersome and not very 
illuminating, so we just give the final results here term-by-term. We abbreviate: 

s^ = sin^^(,, cw = cos9w. (3.105) 

C.VVV --^{\ {W-W+>'W-W+'' - W-W->'WJ-W+'') (3.106) 

+ cw {W^W+^Z^Z'' - W-Z'^W+Z") 
+ {W-W+^'A^A" - W-A^W+A") 
+swcw {2W-W+^'A,Z'' - W-Ai'W+Z'' - W-Z^'W+A'')} . 

= -I— [((^1, - d2p) g^" + {d2, - ds^) g^P + (^3. - 5i.) gn (3.107) 

Sw 

xW+''{xi)W-^{x2){cos9^Z''{xs)+sm9^A'^{xs))] 
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the coordinates {xi, x^, X3) are written out explicitly here only to notate that the deriva- 
tives ^2, ^3) in parentheses act on different fields. 



e2 



^vvss -— {:iA,W+''hx- - iA^W-^hx-^ - A^W+^xsX~ - A^W'^x^X^) (3-108) 



2sw 

„2 



+ {-iZ^W+^hx- + iZ^W-'^hx^ + z^w^i'x^x- + z^w-i'x^x^) 



+ e'A.A^x-'x- + "^^^^^A.Z'^x^x- + "^^^Z.Z^x^x 



+ ^ {w;w-^hh + w;w-^x^x^ + 2w;w-^x'x-) 

( 



+ {Z.Z'^hh + Z.Z^x^x^) ■ 



C^vs ^emw {iW^A^x- - iW' A^x-") + ^^^^^ {-iW;: Z^x' + iW;: Z^x") 

+ "^W^W-^h + ^^Z.Z^h (3.109) 

Sw 2swCyy 



^vss -— [W+'^ihd^X- - X'd^h) + W-^ihd^x' - x+a^/i)] (3.110) 



2sw 

ie 



2s 



w 



[W+'ix^d.X- - X'd.Xs) - W-^{x3d,x^ - x+^^Xs)] 



+ieA^(x~d,x'' - X^d.x-) + ie^^-^Z^ix-d,x^ - x^d.x') 
+ ^^^^'^(/i5^X3-X35^/i)- 



2 2 

^ssss = - 00 2 "2 {hhhh + X3X3X3X3 + 4x^X^X"X~ (3.111) 



^2syyTnyy 



+ 2hhx3X3 + 4:hhx^X +^X^X X3X3) 



^ (/i/i/i + /iX3X3 + 2/iX+X") ■ (3.112) 



Note that we expressed the vacuum expectation value v and the constant of the Higgs 
potential A in terms of the physical masses: 

2mH^sin6'^„ /i^ [ niHe V /oiio\ 

V = , A = ^ = — 1= . (3.113) 

e \2V2mwsme^J 
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3.3.2 Cgfi - non-linear gauge-fixing 

We now switch to the analysis of the interaction terms arising from the non-linear gauge- 
fixing Cpfi. As we discussed when analyzing the gauge-fixing in general, the constraint 



functions can be of the form (2.34) 



G'' = ^,A''>'-K\x^ + G':,,, (3.114) 

with G^j being any quadratic functions of the fields, and the gauge-fixing Lagrangian is 

^gf = -^G'U"^G'^ (3.115) 

The quadratic additions G^j obviously give then the interaction terms. We have actually 
already solved for matrices L and K that give nice propagators in earlier sections, and so 
we can rewrite £gf after changing to the basis of physical fields as (keeping C^i undefined) : 

£gf = - ^G-G^ - Ug'? - Ug"'? (3.116) 
-^{d,Z^-^,mzx' + Gl,f 

1 /„ ... ^^x2 



{d.A^ + GtiY 



2U 

where the non-linear additions were transformed into G^j (and its hermitian conjugate 
G~j), which can be a complex quadratic function of the fields, and Gfj, G^j that are real 
quadratic modifications. It has to be kept in mind when choosing the quadratic additions, 
that, in order not to introduce unnecessary complications, the resulting £gf should respect 
the unbroken symmetries of the Lagrangian, which are the electromagnetic U{1) and the 
5'f/(3) of QCD. That means that G^j should transform as a positively charged field, Gfj 
and G^i should be neutral, and they all should, of course, be neutral under SU(?>). 

In this work we will assume a particular class of non-linear gauges given by the choice, 
in accordance to [5], [S]: 

G+ = ~^&e{A,W^^) - ^/3^(Z,W^+) - 5^(/.x+) + m^(X3X+) (3.117) 

Gl, = -e-^{hx,) (3.118) 
Gt, = 0. (3.119) 
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We have introduced here five new gauge-fixing parameters 

a, (3, S, K, e 

that specify the concrete choice of the gauge. The case when they all are equal to 
corresponds, of course, to the usual linear gauge, and their variation allows to explore the 
different non-linear gauges, and check the results of various calculations for more general 
gauge-invariance . 

The remaining goal of this section is now just to plug in the definitions for G^j and 



Gfi into (3.116) and to collect the terms of higher than quadratic order, which will give us 
£gfi. Remember that the quadratic terms in £gf are independent of G^j and were already 
taken into account for calculation of the propagators. Again, we split £gfi into pieces 
describing the interactions of different particle species: 

£gfi = C 

eivvvv givvv ~t~ ^givvss ~t~ ^givvs ~t~ ^gfvss ^gfssss ~t~ ^gisss^ (3.120) 

and give the results term-by-term. There are no peculiarities involved in this calculation 



-just the straightforward squaring of parentheses in (3.116): 



^gf.... = - -^{W-A^W;^An - ^—^{W-Z^W^Z'^) (3.121) 



-iy^-Am^z'' + w-zm^A"") 



Hgivvv ~{W^A^d,W~'' - W-Af^d^W^n + "p^iW^Z^d^W-" - W-Z'^d^W^"), 

(3.122) 



C,ivvss = - {A.W^^hx- - A,W-^hx^) + [A.W^^x^x- + A^H^-^X3X+) 

_ 'l^^z.w^^hx- - z,w-^hx^) + ^^f^{Z,W-^^X3X- + Z,W-^X3X' 



£gf,,, = - laemwiW;: A^x- - ^-A^x^) - 'J^^^^I^^iw; Z^x' ~ W'^Z^X^). 



(3.123) 



(3.124) 
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C-^ivss =^{hx-d,W+'^ + hx^d,W~n + 7^{X3X-d,W^^ - Xzx'-d.W-'^) (3.125) 



+ 



ee 



glssss 



4,2 



-{hhx^X ) - 



4,2 



iX3X3X'^X ) 



fV6 

8s2 p2 



ihhxsXa) 



Cgfsss = : [hx X ) - ^„ „ [hXsXs)- 



(3.126) 



(3.127) 



All these gauge-fixing terms combine with the regular terms from £; to effectively alter 
the strength or the structure of the coupling, depending on the gauge parameters. We 
will show that when we write out the full Lagrangian with all its terms. 



3.3.3 £ghi - ghost interactions 



We now move on to the last part of the interaction Lagrangian £ghi - the ghost interaction 
terms. Remember from the description of the gauge-fixing procedure that the ghost 



Lagrangian is given by (2.19): 



gh 



(fix) 



SG'^lA^.Xc 



6a^ 



0=0 



c\x), 



(3.128) 



and that we have already diagonalized the ghost propagators to get the ghost energy 



eigenstates as (3.66) 



V2 



(cr\ 



CwC^ — SwC^ 
\SwC^ + CwC^J 

/^(ci + .c2)\ 
^(c--.c-) 

CwC^ — SwC^ 



(3.129) 



(3.130) 



Note also that we had to multiply the variation 6G/6a by (—g) to get (2.31 ) and correctly 
normalized ghost propagators g We can then express the ghost Lagrangian directly in 
terms of the physical fields as: 



6G^ 
6a^ 



db^rb 



-c"" (^S) (3.131) 



^It should be kept in mind that (g) is not really a single constant, but a diagonal matrix g"^^ with 
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with transformed and 



a 



(3.132) 
(3.133) 



The transformed constraints G'"" are, of course, just the physical constraints as in (3.116): 

/ gA I d,W+^ - ^wmwx'' + G+A 
G- d.W-f^ - ^wmwX- + 



G" 



G^ 



\ 



(3.134) 



and the transformed a'" are the physical gauge transformation coefficients in terms of 
which we will have to find variations of the fields. So, actually, all we have to do now is 
to calculate the variations of the physical fields in terms of a'"", and then we can simply 



read off the ghost Lagrangian using (3.131 ) and the definitions for G , G , G . 



The variation in the physical gauge fields is: 



(3.135) 



After transforming the structure constants we get: 



I" \ I" Sw 

6Z^ =d.a^ - ie— (W+a~ - W-a+) 



cw 

Sw 



6Ai, =df,a' 



a 



(3.136) 
(3.137) 
(3.138) 



Note that the new gauge transformation coefficients are indeed physical in 

the sense that each of them yields a transformation of the corresponding physical gauge 
field by df^a. The non-Abelian part of the gauge transformation naturally doesn't look 
symmetric anymore when written out in terms of the physical fields, because the symmetry 
is actually broken. 

Now, to get the variations of the scalar fields we similarly need to transform everything 
to the physical basis. Here it's convenient to use a notation similar to that for gauge 
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bosons, combining Goldstones and Higgs into one vector: 



X~ 

X3 



(3.139) 



with (looking back at (3.59)): 





1 

V2 


i 








\ 




1 
V2 


i 


















1 







1 











1 


/ 



(3.140) 



The variation is then: 



which results in 



5X3 

5h 



mzce 



2sw 
ie 

2sw 



X^tt" - X"«" 

o , X^tt" + X^a^ + — Xstt^ 
zsw V '^w 



Cw 



(3.141) 

(3.142) 
(3.143) 
(3.144) 



Note that again the terms just proportional to a end up being diagonal in the physical 
basis, which is just a check that we did the transformations correctly: their diagonality 
is equivalent to the ghost propagators being diagonal, and that's what we started with in 
this section. 

We can now use these expressions for variations to write out the ghost Lagrangian. 
First consider the part of Cghi which is present even without the non-linear gauge-fixing 



terms. Plugging in (3.134) into (3.131) with G^i = yields: 



ghO 



ba 



la 



(3.145) 



ba 



ba'' 



5X3 

ba"" 



Note that the operator 



ba'' 
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just replaces the a'" in the variations with c'". We can then easily write out (3.145) by just 
substituting the variations of all the fields. Also, since we are looking for the interaction 
terms, we drop the quadratic ghost-antighost terms, yielding: 

.Cghio = + ^e{^^c-) ( W+^c^ + ^W+^c^ - A^c^ - ^Z^cA (3.146) 

\ Sw Sw J 

- ie{d.c+) ( W-^c^ + "^W-^c^ - A^c- - ^Z^c'] 

\ Sw Sw J 

Sw 

+ c- ('i^^^c^ + 2swX^c^ - + ^hc^ 

^Sw \ Cw 

c+ X c + 2swX c -xsc - the 

■^Sw V Cw 

ie^zmz -z f + - - + ^ . ^ 
c X c — X c he 



2sw \ Cw 

This Lagrangian describes the ghost-ghost-vector interactions (first 4 lines) and ghost- 
ghost-scalar interactions (last 3 lines) and it would be the final result in the case of linear 
gauge-fixing. 

With non-linear gauge-fixing the ghost interaction Lagrangian will also contain terms 
arising from the non-zero G^j, Gfj: 



These terms can also be calculated straightforwardly using the definitions (3.117). It's 
convenient to work out the variation of the constraints term-by-term, which gives the 
Lagrangian split up in terms, each proportional to a different gauge-fixing parameter: 

•^ghinl -^ghio ~l~ ^ghif3 ~l~ ^ghiS ~l~ -^ghi/t ~l~ -^ghif:- (3.148) 

We list here the final resulting expressions: 

Cgua = + iaec~ (A^9^c+ + W+d^'c'^) - ae^c' (^A^W+f'c^ + ^A^W+^'c^- 

-A^A^'c^ - ^A^Z^'c^ - W+W+^'c- + W+W-^c+^ + h.c. (3.149) 
Cgkif^ = + i(3e^ c- {Z^d^'c+ + W+d^c^) - /^e^f^ c" (^Z^W+f'c^ + ^Zf.W+^'c^- 

-Z.A^'c^ - '^Z^Z^c^ - ^W+W+^'c- + ^W+W'^c+) + h.c. (3.150) 
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C,M5 = - S^^^-^ c-hc^ + 6-^ c- (x+X+c- + X+X-C++ (3.151) 
Zsw 4:S^r 

1 p2 _ 2 \ 

+— X^Xsc^ + i— -X^hc^ + 2iswX^hc^ - ixshc^ - hhc+ + h.c. 



Cw Cw 

>Cghi« ^ii^T^ c" imwX^c' + mzx^c^) + '^-y^ c" f-x^X'^c" + x"^X~c++ (3.152) 
+—X^hc^ + ^^~^^ ^+^3c^ + 2sw^x^X3C^ + ix^hc-^ - X^X^cA + h.c. 

Cw Cw J 

„ -,e^zmz_zu z , ~ ^"^^z _z /.+,_,._,+ , + _ , _ + , 

>Cghi£ = - e- c he + £— ^ c [-IX he + ix hc^ + x Xac + x X3C^+ 

2swcw 4:sf^cw 

+— X3X3C^ - — /i/ic^) . (3.153) 

Cw Cw J 

The abbreviation "h.c." above means that a hermitian conjugate of the whole expression 
has to be added as well. Among the terms above we see some ghost-ghost- vector and ghost- 
ghost-scalar interactions, which provide modifications to the terms in £ghiOi but also we 
have lots of new terms describing ghost-ghost-vector-vector and ghost-ghost-scalar-scalar 
interactions that are only present in the case of the non-linear gauge-fixing. 

For the final ghost interaction Lagrangian we just need to collect all the parts fisted 
in this section: 

■^ghi = ^ghiO + -^ghia + -^ghi/? + ^ghii5 + -^ghiK + ^ghi£. (3.154) 



3.4 Full Lagrangian 

In the previous sections we have worked out all the terms needed to construct the full 
Lagrangian of the Standard Model in the case of non-linear gauge-fixing. We provide here 
a summary of the contents of the Lagrangian and it is written out explicitly in Appendix 
A. The full list of interaction vertices arising from the Lagrangian is given in Appendix 
B. 

The Lagrangian contains the following pieces: 

+ + + (3.155) 

WW wv wss vvs vss ssss ~l~ ^sss'^ 
~l~ ^ggv ^ggs ^ggvv ^ggss- 

The terms C^^CsiCf^Cg hold the quadratic terms for vector, scalar, fermion and ghost 
fields respectively, describing their propagators - fermion propagator terms can be found 
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in Section 3.2 and the others in Section 3.1 The rest of the Lagrangian describes the 
interactions: the terms not involving fermions or ghosts are found by combining the results 
of Sections 3.3.1 and 3.3.2, the fermion interactions £//^,£//s are described in Section 



3.2| and, finally, the ghost interactions were listed in the previous Section |3.3.3 



We have now derived the full Standard Model Lagrangian, which we will apply now in 
calculations. First, in the following section, we will give an example of a tree- level process 
amplitude calculation, which explicitly demonstrates the use of the derived Lagrangian 
and the effects of the non-linear gauge-fixing. Then, finally, we will describe the imple- 
mentation of the Lagrangian in FeynArts and the calculations we did using the package, 
which is the real purpose of this work. 



4 Explicit tree-level 77 WW calculation 

In this section we will perform a tree-level calculation for the tree-level process of W^W~ 
production in a photon-photon interaction, using the non-linear gauge-fixing Lagrangian. 
The goal of this calculation is to see, that all the gauge-fixing parameters actually cancel 
out in the final result as expected. This particular process was chosen as an example 
because it illustrates well the non-trivial cancellation of the gauge-fixing parameters. 
The tree-level diagrams for this process are the following: 



39 




and the full amplitude is: 

M=Mu + Miu + M2t + M2u + M3. (4.1) 

The interesting fact we will see is that each diagram separately is not independent of the 
gauge-fixing parameters, corresponding to the fact that each diagram separately is not 
gauge-invariant, and only the sum of all 5 diagrams gives a gauge-invariant (and thus 
independent of the parameters) result. 

We now proceed to the calculation of the diagrams - all the needed Feynman rules for 
our Lagrangian can be found in the Appendix B. Note that all the interaction vertices 
involved in the process depend on the non-linear gauge-fixing parameter a and our interest 
will be how it cancels out in the final amplitude. 

Let's define: 

Qt =Pi-P3=P4-P2, t = ql = - 2pi ■ = - 2p2 ■ Pa, (4.2) 
Qu=P4-Pi=P2- P3 > u = ql = mw-2pi-p4 = mw- 2p2 ■ ps, (4.3) 
where t, u are Mandelstam variables. Also note that 

Pl-P3=P2-P4, Pl-P4=P2-P3, (4.4) 

and clearly 

= P2 = 0, pI^pI^ m^. (4.5) 

We will use for the polarization of the external vector particle with momentum pi, and 
it is - as usual - orthogonal to the momentum of the particle (no sum over i here) : 

Si^-pf^O. (4.6) 

Throughout the calculation the following abbreviations will be used: 

S'ij = S'l * or 

eij=ei-ej, or (4.7) 
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there is no potential confusion here, because each polarization vector appears either as £j 
for incoming particles, or as e* for outgoing, but not both. 

One final note on how the calculation will be organized. In order to analyze the 
dependence of the final result on a it is convenient to separate terms in the amplitude 
containing different orders of a. There will be terms of first and second order in a, which 
we will notate as: 

M =Mo + Ma + Ma2, (4.8) 

Mo ^Mito + MiuQ + M2tQ + M2U0 + MsQ, (4.9) 

Ma ^Mlta + Mlua + M2ta + M2na + M^a, (4.10) 
Ma2 =M + M ^M 2ta2 + M2ua2 + M^a2- (^-H) 

We will proceed by calculating each diagram in a row and splitting the result into three 

terms by the order of a for each diagram. 
The first diagram is now: 

Mit = - ieixe2aelyel^C^^^{-qt, -P3,Pi)Cwwa(-P4, quP2) 



X 9,, - (1 - iw) ^ ■ (4-12) 

We use a shorthand Cwwa(Pi>P2,P3) for the WW A interaction vertex (see Appendix B), 
where (//, A) are contracted with incoming W~ , W'^ , A lines respectively, and (pi,P2,P3) 
are their respective momenta. The interaction vertex Cwwa consists of two terms: Cqwwa, 
the part coming from the original Lagrangian, and the gauge-fixing part C^wwa which is 
proportional to a. Also note that 

Cwwa(Pi>P2,P3) = -C^wa(P2,Pi,P3)- (4.13) 
Let's rewrite Mu then as: 

Mu -^1A£2.£3.^4.^37;^ ^9,P - (1 - ^W)^_^^^2^) 

X (CowwA(-^t' -P3,Pl) + CgWWA(-?t' -P3,Pi)) (4-14) 

X (CoWWA(5t' -P4,P2) + Cg^AiQt, -P4,P2)), 

and start with the term quadratic in a, which will come from taking Cg for both interaction 



41 



vertices. It is: 



X 









)■ 








ea 


■( 






t 



3 • ' vv y 



2 



(4.15) 



Here the third equahty is from e ■ p = 0, the fourth one from = t, and the fifth one is 
just some simphfying algebra. 

Now we move to the term Aiua proportional to a, which comes from (4.14) when 
choosing Cg for one interaction vertex and Cq for the other: 



X (CowWAi-Qt, -P3,Pl)C^wWA(9t' -P4,P2) (4.16) 
+ C'oWWA (^t ' , P2 ) CgWWA ( > -P3 , Pi ) ) • 

Vertex Cg again reduces to just a constant times qt because of orthogonal polarizations: 

(1^) - (1 - 6.)^^) (4.17) 
x(CowwA(-9t' -P3,Pi)gtV - C'owwaI^*' -P4,P2)qt9''^)- 
In order to evaluate this expression we need to calculate the piece: 

-^13 =(ltfj.£lx£*3uCowWAi'~1t^ -P3,Pl) 

= {Pl - P3)m^1a4^C'oWWa(P3 - Pi, -P3,Pl) 

= - 2e(pi - P3)m^ia£L ((-2P3 + Pi) V" + (Pi + P3)^^7'^' + (-2pi + PsTg^'') 

= _ ,e {-2{p^ ■ sDips ■ e^) - {pl-p\){e, ■ el) + 2{p, ■ e^){p, ■ el)) 

=iem'^ei3. (4.18) 
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There is also an analogous contraction: 



K24 = — (ltti^2\^\u^QWWpS.<it, —Pi,P2) 

= {P2 - P4)m£^2A£L<^0WWa(P4 - P2, "^4,^2) 



which, comparing to the second line of (4.18), is clearly the same as K13 but with all 



indexes 1 replaced by 2 and 3 replaced by 4. Therefore, the same calculation yields: 

^24 = ^em^e24- (4.19) 
Then Aiua is just repetitive occurrence of Ki^ and K24: 



. , tea 

Mlta 



ieoi 



Kl^e2i + i^24£l3 - ^ (-^13^24^4 + K24,ei3qt) 



iw{t - m 



2 2iemy^ei3e2i 



w! 



t - iwm^w. 



m 



= - 2e'ae,se24 , , 2 ^ (4-20) 
t — c^w'mw 

where the last line is again just some algebra. 

Finally, the last term in A^ii, independent of a, is with both vertices taken as Cq: 

Mm =eiAe2a£L^LCowwA(-^t' -P3,Pi)C^wwA(?t> -^4,^2) (4.21) 



We would like to do one last separation and consider the piece of this amplitude: 

-^ItOO =^lX^2a^*3u^lKC^0WWAi'~^t' ~P3) Pl)C'oWWA('?i; ~Pi^P2) , ^^''2 ■ 
f ■ \2 * * 9 up 

=(-^e) eixe2.e^^e^ 

X ((2p3 - Pi)^9'"' - iP3 + PiTg"' + (2pi - PsYg^^) 

X ((2p4 - P2r9'"' - iP4 + P2Yg'"' + (2p2 - PAYg"') . (4.22) 

The reason to write it separately is that this piece is completely independent of the gauge 
parameters (unlike the rest of M-uq-, which has ^w)- This is in fact the full M.it in the 
case of simple gauge choice = 1; = 0, and we will see in the end that it contains 
most of the full amplitude in the t-channel. However, in the context of this work this term 
is not very interesting, because we want to analyze the dependence on the gauge, and it is 



already independent. Furthermore, expanding (4.22) into scalar products is quite tedious 
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and not very illuminating, so there is not much more we can do, but to keep this term as 



A^itoo- The rest of A4uo reduces to Ki^ and and we have from (4.21): 



We have now calculated the full M-u-, which can be built from (4.15), (4.20) and (4.23). 
The next diagram on the list is M.\u- Here, however, almost no extra calculations are 
needed, because this amplitude is exactly like M.\t but with the two photons interchanged. 
That is, we can rewrite completely the same calculation making everywhere the following 
replacements: 

Pi, (4.24) 
£1, (4.25) 
I. (4.26) 

Therefore, the final results are analogously: 



V\ - 


-^P2, 


V2 


e\ - 


^2, 


£2 


It - 




t - 



u 



Mxuoa =e a 61^62^1—, r (4.27) 

.2 



and 



Miua = - 2e^aeue23 r^^, (4.28) 

Miuo =-Mi«oo - 6^6146237 ""Hl'^r^ 2 ^ ^ (4-29) 



ki /•\2 ** 9 HP 



u - 

X ((2p3 - P2)'9'"' - {P3 + P2r9'' + {2P2 - P3r9'') 

X ((2p4 - Pir9'^ - {PA + P1Y9'"' + (2pi - PAr9'") . (4.30) 

Now we move to calculation of M.2t and M.2u- The interaction vertex with a Goldstone 
boson and its propagator are quite simple, so we can write the amplitude out fully as: 

M2t = - iei^e2pel^el^{-{l - a)emw9'"'){{^ - a)emw9'"')-. — ^ r 

t — i^w^w 

2 

= - e'(l -2a + a2)£i3£24-— (4-31) 

t — i;wmw 
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The ti-channel part is analogously 



,2 



M2U = -e\l -2a + a^)eue23 r^^, (4.32) 

u — l;wn^w 



and we can write out components with different powers of a: 

2 2 

M2to = -e^£i3£24T — 2 ' -^suo = -e^£i4£23 2 (4-33) 

1 1 
TTX 771/ 

M2ta = 2e'a£i3£24 , 2 ^ -^S^a = 2e'a£i4£23 (4.34) 

2 2 

M2ta2 = -e^Q;^£i3£247 2 ' •^2«a2 = -e^Q;^£i4£23 ^ 2 (4-35) 



And, finally, the last diagram AI3 is simply the interaction vertex with two VF's and 
two photons: 



Ms = - iei^e2u£lp£lJe^ 



a" 



— 2£i2£34 + I 1 ~ 7 — (£13^24 + £14^23) 



(4.36) 



The (5-independent component is 

Mm = (-2£i2£34 + £:i3£^24 + 14^23) , (4.37) 

and it is convenient to further split up the part quadratic in a into two terms - corre- 
sponding to i-channel and li-channel respectively: 

Mzta2 — T £13^24, Mzua2 — T £l4£23- (4.38) 

The full M^ is then 

Mz = M30 + MMa2 + M^ua2, (4.39) 

and there is no term in M.^ proportional to ct, that is, M'ia = 0. 

We have now all the parts to construct the full amplitude M. We will start with 
a-dependent terms and show that they add up to 0. Then we will analyze the Mq 
component and show that the other gauge parameter, ^w, also drops out. Since we have 
split Al3a2 into t and u components, we can perform a-dependence calculation for t and 
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u channels independently. So the quadratic t component is: 

Mta2 + M2ta2 + Mua2 

2-2 ^ 2-2 e^a^ 

— e a £13^24 7 — 77 ^ 2^ — e a £l3£24 7 ^ 2 1 £13^24 

_g2^2^ e ( - ^ 

2-2^ o ^ ~ ^wrri^ - t + iwm\r \ „ / , 

" i iw(t-i^^i,) ) = °' 

and by complete analogy 

Mua2 = Mlua2 + A^2«a2 + A^3ua2 = 0, (4.41) 

the quadratic a terms have canceled! The linear terms also add up to as expected: 

2 2 

Mta = Mita + M2ta = -2e^a£i3£24 , 2 + 2e'fei3£24 , 2 = 0' (4-42) 

and so also 

Mua = Miua + M2ua = 0. (4.43) 

So the full amplitude is a independent! It is explicit here how only the sum of all diagrams 
gives this gauge invariance, while each separate diagram is gauge dependent. 

Now let's finish with calculating A^o, where we still have to see the {vF-independence. 
Since already has no we evaluate the sum of the other two terms: 



MuQ + M2m =Alitoo - e £13^^247^ ttt, 7 ^ - e £i3£24 



{t - m^) {t - iwm\^) t - ^wrn^w 



^2 /rn\r{l-iv 



— Alitoo — e £13^247 — \ — r i 



t-Cwm^ \ t-m 



--Mitoo - e^£i3£24T 'V' (4-44) 



and likewise for w-channel: 



Miuo + M2U0 = A^i„oo - e^£i4£23 (4.45) 



The dependence on has dropped out! The full amphtude is then 

M =Muo + M2to + Miuo + M2U0 + M30 

2 2 
2^ ^ '^w 2^ ^ "^w 



— ^UOO + -^luOO ~ e €13624- n e £i4£23 



t - u - 



+ (-2£i2£34 + £13^24 + £l4£23) , (4.46) 
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which is free of any gauge-fixing parameters. As mentioned before, the largest part of 



the result is contained in TWitoo and TWiuoo defined in (4.22), (4.30), and for our purposes 
there is no need to further expand them. 

In this section we have performed a full tree-level calculation of the amplitude for 
W~^W~ production, and showed that it is independent of both non-linear and linear 
gauge-fixing parameters. On the other hand, it is clear now how this fact can be used 
to simplify the calculations. Choosing, for example, the gauge a = = ^ ehminates 
the diagrams A^2t and A42u, and removes t and u channel components from Ai^ - the 
only term remaining there is proportional to 612634. Then the diagram Aiu with the 
modified interaction vertex will produce the full amplitude for the t channel. Specifically, 
in addition to Aiuoo, which is the "usual" amplitude for that diagram, corresponding to 



gauge = I, a = 0, it will produce the other two terms in (4.46) proportional to 61-3624 
that are normally seen as coming from the diagrams A^2t and Ai^. 



5 Non-linear gauge-fixing in FeynArts 

We have implemented the full Lagrangian with non-linear gauge-fixing in FeynArts, and 
used it to perform a number of automatic amplitude calculations. The results, as expected, 
turned out to be gauge-independent, which is a good indication that there were no errors 
in the implementation. 

In this section first we will briefly describe FeynArts/FormCalc and how it is used. 
Then we will explain what changes have to be made to accommodate non-linear gauge- 
fixing, and, finally, we will discuss the calculations that we did. 



5.1 About FeynArts and FormCalc 

FeynArts/FormCalc P] is a freely available open source package for Mathematica that 
performs automatic calculation of Feynman diagrams up to 2-loop level. The download 
of the package and more information about it is available at: 

http : // www . f eynart s . de 

In this work we have worked with the version FeynArts 3.2. 

The calculation of an amplitude in FeynArts is performed as follows. The process is 
defined simply as a set of incoming and outgoing particles. The package then provides 
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functions to do all the required steps: first, it generates all the possible Feynman diagrams 
at the given loop level, including the counterterm diagrams, second, it constructs the alge- 
braic expressions for the diagrams using the definitions for vertices and propagators, and, 
finally, it can algebraically evaluate these expressions by applying various rules resulting 
in a greatly simplified symbolic result for the total amplitude. During this last step the 
loop integrals are reduced to the Passarino-Veltman integrals, also the counterterms are 
evaluated using the defined renormalization scheme, so the resulting amplitude should be 
already free of ultraviolet divergences. 

The final symbolic expression for the amplitude can be used then to calculate numerical 
results. For that purpose the package can generate a full FORTRAN program, which, 
when compiled and run, will loop through the phase space of a process producing numeric 
values for the cross-section, which can be compared with experimental results. In this 
work, however, we have not used numeric calculations. 

The first steps - generating diagrams, and translating them into expressions using 
propagators and vertices - are done using the functionality of Mathematica and this 
part of the package is called FeynArts, while the last steps - algebraic evaluation and 
FORTRAN code generation - internally uses another program FORM, and so this part 
is called the FormCalc package. Nevertheless, these two are closely related and can be 
really thought of as one complete package. 

The really attractive feature of FeynArts is that the physical model for calculations is 
not built-in, but provided by external input files called model files, which contain all the 
information about the theory: the definitions of the fields, their properties, propagators, 
interaction vertices, counterterms and renormalization conditions. The package comes 
already with some predefined models, including QED, Standard Model and MSSM, that 
can be used immediately for calculations. In order to implement the non-linear gauge- 
fixing we had to modify the existing Standard Model file, and we give now the details of 
the needed modifications. 

5.2 Implementation in FeynArts 

The existing definition of the Standard Model in FeynArts matches the one we have 
discussed in this work in the case of linear gauge-fixing. That is, all the fields, propagators 
and interaction vertices are the same as listed in Appendix |b] with the choice a = (3 = 6 = 
K = e = 0. Therefore, the only change needed for non-linear gauge-fixing is to modify the 
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Table 1: Extra factors to get a FeynArts vertex from the one in Appendix [B] 



Field 


Factor 


Z 


-1 

+i 
—i 
-1 

respectively 



interaction vertices with the additional terms proportional to the non-linear gauge-fixing 
parameters. 

First, in order to compare our vertices with the ones defined in FeynArts, we note 
the following differences in sign conventions. FeynArts defines the Fourier transform so 
that 9^ ipfj^ as opposed to (9^ —iPtJ., so there is a relative minus sign in all vertices 
that include particle momenta. Other differences relate to the additional factors in the 
definitions of the fields (which are, of course, allowed), as a result of that vertices get 
multiplied by a factor for each occurrence of the fields with differing definitions. We give 
in Table 1 these extra factors that should be applied to a vertex given in Appendix [B} 
in order to get the corresponding FeynArts vertex. In order not to be forced to change 
all the existing definitions for interaction vertices in FeynArts, we adopted the FeynArts 
conventions for the fields when implementing the non-linear gauge-fixing modifications, 
thus converting the interaction vertices using these rules. 

Consider now the new terms in the interaction vertices that the non-linear gauge-fixing 
yields. Note that some of them just change the overall factor of the already existing 
vertices, like for Vector- Vector-Scalar-Scalar or Vector- Vector-Scalar couplings, while in 
other cases the non-linear gauge-fixing gives a whole new kinematic term, e.g. for Vector- 
Vector- Vector or Vector-Scalar-Scalar couplings. The latter case is especially notable in 
ghost interactions, where you can see, that the whole two new classes of interactions - 
Ghost-Ghost- Vector- Vector and Ghost-Ghost-Scalar-Scalar - were generated by the non- 
linear gauge-fixing. Now the reason it is important to note the new kinematic terms is 
that they require some additional work in order to be implemented. That comes from 
the fact that FeynArts defines a physical model in two steps: the generic model file - 
called Lorentz . gen in the case of the Standard Model - and the classes model file - 
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called SM.mod. The generic model file does not refer to the specific particles, but only to 
generic fields, which can be vector, scalar, fermion or ghost. The file then defines generic 
couplings saying that, for example, the Vector- Vector- Vector vertex is of the form: 

CryV(Pi>P2,P3) = C'yyy ■ [g^'' {p, - p^Y + 9"" {p2 - PsT + 9''{P3 - , (5.1) 
and Vector- Vector- Vector- Vector is: 

^vvvv ~ ^vvvv y y ^ ^vvvv y y ^ ^ww y y i 

so that just the coefficients C^, C^, depend on the fields, but not the structure of the 
coupling. The classes model file then defines exactly what kinds of fields there are in the 
theory, and specifies these coefficients for all specific vertices. 

This separation into generic and classes models is, of course, very useful because it 
allows to define common properties for the vectors, scalars, fermions, ghosts, and thus 
the same Lorentz.gen file can be used, for example, for both SM and MSSM. The 
complication with the additional kinematic terms arising from the non-linear gauge-fixing 
is that they do not "fit" into the standard Lorentz . gen definitions. For example, the 
generic Vector- Vector- Vector vertex has to be defined as 

C^v'v{Vi.P2,Pz) =Cl.yy ■ [y^'ipi - P2Y + y'''{P2 - PzY + y"' {Pz " Pif] (5-3) 

+ Clyy ■ p^.g^" + C^yy ■ p^,g^^ + C'yyV " 

in order to accommodate the extra terms proportional to Piy"^^ — pl^gPP. Note that Cyyy 
is not really used in the vertices that we have, but it has to be included in the generic 
coupling definition, because FeynArts requires that the generic coupling keeps the same 
structure under the permutation of the fields. The other generic couplings we had to 
extend are: 

C'ssviPi^P2,P3) =C' ■ {Pi-P2r + C^-Pl (5.4) 
^ggvv ~ ^ggvv ' 9'^" ^ (5-5) 
Cggss = Cggssi (5-6) 

where in the first one we had to allow for the new extra term proportional to pg, and the 
following two are the four-vertices with ghosts, that were not present at all. The parts of 
Lorentz.gen containing these changes can be found in Appendix [Cj 

Once the changes to the generic model file are made, the modifications of the class 
model file are quite straightforward: we just have to go through all the vertices, and 



50 



include the terms proportional to a, /?, 6, n, e into the existing or newly defined coupling 
constants. Some examples of these alterations are also given in Appendix [C| 

5.3 Calculation tests 

After we have implemented the non-linear gauge-fixing by modifying the generic and 
class model files, we have used the new model to perform calculations of a number of 
amplitudes in order to test the implementation. The condition for testing is, of course, 
the same fact that this whole work relies on: the final amplitude for a process has to be 
gauge-independent, which means that in the calculation of any physical process the newly 
introduced gauge-fixing parameters a,/?, 5, e have to cancel out. Here we describe the 
tree-level and loop-level calculations for which we have checked this cancellation. 

5.3.1 Tree-level calculations 

For the tree-level amplitude test we have chosen to do the same 77 W^W~ calculation 
that we already did by hand in Section |4| in order to be able to analyze it better. The 
immediate result for the amplitude after writing out the diagrams in terms of vertices 
and propagators looked highly dependent on gauge-fixing parameters and a, which 
is to be expected, since, as we saw, all the vertices involved do depend on a. However, 
after running FormCalc, which is able to simplify the result by using the fact that the 
polarizations are transverse, and then using simplifying functions of Mathematica the 



gauge-fixing parameters did cancel out and the final result corresponded to (4.46). 

By using the same simplifying procedures we were also able to confirm the gauge- 
invariance of other tree-level processes including ZZ W~W~^, e~e^ W~W~^, 

5.3.2 One-loop calculations 

The one-loop calculations are substantially more difficult than the tree-level, and we have 
not been able to show algebraically that the full one-loop amplitude is independent of 
the gauge parameters. Note that a single one-loop amplitude can contain almost all the 
vertices of the theory in its internal loop, therefore, the cancellation of all the different 
gauge-fixing terms is far more complicated than in the tree-level (which already was non- 
trivial), and the gauge-invariance check is much more strict. Such a check would likely 
require a numeric calculation using the generated FORTRAN program, which can be run 
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using different choices of the gauge-fixing parameters, and the results should not vary. 
However, before a numerical calculation can be carried out, it should be checked that the 
amplitude is free of ultraviolet divergences, otherwise the result of the numerical calcula- 
tion (which doesn't include the divergent terms) could not be trusted. The cancellation 
of the divergent terms has to be ensured by the renormalization scheme, and for the 
definitions in the given SM . mod, which has been used extensively in one-loop calculations 

it has been tested to work. But now, after modifying the model with the non-linear 
gauge-fixing, we have to first check again that the renormalization works and the UV 
divergences cancel before doing any further one-loop calculations. 

FeynArts provides a way to check for UV divergences using the function UVDivergentPart, 
which returns the coefficient of the divergency of the final amplitude. If the renormal- 
ization is successful, this coefficient should be zero. We have performed this check on 
the one-loop amplitudes calculated with the non-linear gauge-fixing, and, using the sim- 
plifying functions of Mathematica, showed that all the tested amplitudes are indeed free 
of UV divergence. The one-loop amplitudes that were tested include all the two-point 
functions and the processes e~e+ W~W^, e~e~^ f^'f^'^, 77 W^W~ . Thus we 
have shown that the renormalization scheme of the Standard Model works when the non- 
linear gauge-fixing is included, and the numerical results of one-loop amplitudes can be 
considered. 

6 Conclusions 

In this work we have presented a full derivation of the Feynman rules of the Standard 
Model with the non-linear gauge-fixing. We have implemented these rules in the Feyn- 
Arts/FormCalc package by modifying the existing definitions for the Standard Model. 
We have confirmed that the new definitions give gauge invariant results in the tree-level 
amplitudes, and we have also checked that the one-loop calculations are free of ultraviolet 
divergences. 

The new definitions have not been tested yet to give fully gauge invariant results in the 
one-loop level; such a check would probably have to be done by numeric calculations of the 
cross-sections, since the expressions for one-loop amplitudes seem to be too complicated 
to be analyzed algebraically. Such full check of one-loop gauge invariance would be the 
^FeynArts uses dimensional regularization, so the divergence is 2/(4 — D) 
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most reasonable next step in the continuation of this work. 

Another direction in which this work could be extended is the implementation of 
the non-linear gauge-fixing in the Minimal Supersymmetric Standard Model (MSSM). 
The automatic calculations in the MSSM arc as important as in the SM, because the 
predictions of both models have to be checked against the precision measurements, and 
the non-linear gauge invariance check would be a very useful tool. 
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Appendixes 

A Full Lagrangian 

Wc give the full Lagrangian of the Standard Model grouped in terms based on species of 
the fields (vector, scalar, fermion, ghost) involved. 

+ + + (A.l) 

WW ~l~ ^wv ~l~ ^wss ~l~ ^vvs ~l~ ^vss ~l~ ^ssss ~l~ ^sss~^ 
+ ^ffv + 

H" ^ggv ~l~ ^ggs ~l~ ^ggvv ~l~ ^ggss' 



^v = - W;{-g^^d^ + (1 - Ci^ja'^a'^ - m^^''^)Pr+ (A.2) 

- ^^^.(-^'^'^9^ + (1 - iz)d^d^ - mlgnZu - \A^{-g^''d'' + (1 - iA)d^d'')A. 

1 1 

3 

£/ = ^ ei(i^ - mei)ei + + 'iti(i^ - m„i)iii + - mdi)di 

i=l 



Cww = - ^ {| {w;w+>'w;w+'' - w;w-^w+w+^) (A.3) 

+ (W-W+^Z^Z-^ - (1 - ;3Vev^)PF-Z'^VF+Z^) 
+ (ly-iy+M.A'^ - (1 - aVeH^)iy-AW+A'^) 

+SWCW (2W-W-^^^A,Z^ - (1 - + W^Z^VT+A'^)) } . 

= - ie^ [Z^ {W'd^W; - W+>'d,W;) + W^"" {Z^d^W- - W'^d^Z^) (A.4) 

^W"" {W^f'd^A^, - Af^d^W^) - a/iw{W'^A>'d^W-'' - W- A^'^^W+'')\ 
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C.VSS =^ (^(l - a~5){A^W+^hx- - A^W->^hx^) (A.5) 

-(1 - aK)iA,W+^X3X~ + A^,W-^'X3X'')) 

+ ^ + P6cl^/sl,){Z,W^'^hx- - Z.W-'^hx^) 

+(1 + Pkcl,/s'^){Z,W+''xzX- + Z.W-^xsx"-)) 

+ e'A,A^x^x- + "^^^^^^A.Z^x^x- + "^^f^^Z.Z^x^x- 

+ {w;w-^hh + w^w-^'xzxz + 2Vf;vf- Vx") 

+ ^^^7^ (^^.^'^/i/i + Z.Z'^xzXz) . 



C,,s =emw{l - «) {iW^A'^x" - iW^A'^x'^) (A.6) 

^'^w^w-^h + ^z.z^h 

W+^'ihd^X' - X~d^h) + Sikx-d^W+i') (A.7) 

+w-''{hd^x'^ - x+a^/i) + ~d{hx^d^w->^) 



2s 



w 



1 f 

+ [W+^{X3d,x- - X-d.Xs) + ^iX3X-d,W+n 

ZSw 

-W-'^iXsd^X'^ - x'-d.xs) - ~<X3X^d,W-n] 
+ieA^(x-d,x^ - x'-d.x-) + ie^^Z'^(x-9^X+ - X+^^X") 
[Z'ihd^Xs - XM + e{hx3d,Z^)] ■ 



Cssss = - „^ 2 ^2 (^^^^ + X3X3X3X3 + 4x^X^X"X" (A.8) 

+2(1 + 2e'izml/m]j)hhx-iX3 + 4(1 + 252^iym2^/m|)/i/ix+X~ 
+4(1 + 2/«2{p^m^/m^)x+X"X3X3) 



2 

A.. = - ^JJI^^ (/i/i/i + (1 + 2£ezm|/m^)/iX3X3 + 2(1 + 2~5iww?w/m]i)hx^X~) ■ 

(A.9) 
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3 J- 

SwCw 



+ sin^ OyjCj'^e - | sin^ Oy^wy^u + | sin^ 6't(,d7''d] 



i=l 



= V <i ^= X^i'frie^LeR + uidULdR - rriuURdL) (A.H) 



Z6 — — 

+ ^= X^i-meeRi^L + rUudLUR - madRUL) 

Xsirnee-f^e + nidd'y^d-muU'y^u) 



2sw'mw 
e 

2swn^w 



rueee + ruddd + m„titi) 



i:^^^ = + leW^^" i{d^c-)c^ + ac-d^c^) + ^e^iy+^ f (5mC-)c^ + /3c~a^c^) (A.12) 
+ ieA^" ((9^c+)c- - (5c+9^c-) + ie—Z^" ({d^c+)c- - Pc+d^c-) 

Sw V / 



Sw 



C,,s = + ^^1^ c- f ~ + ^ x+c^ + 2.H^X+c^ - (1 - /^)X3C+ + i(l + ^)/.c+) 



2s\Y V Cw 



-X c + 2swX c - (1 - ii)x3C - i{l + S)hc 



ieizmz _z f _ _ + ~^^^z , _ (^.13) 



2siy V c-w 



-«^^^'^c+ - ~P'^Z^Z^^c+ - (a + /3)^^^Z'^c+- 

-(a + P^)W+W+''c-- + (a + /3 J)iy+iy-^c+) + h.c. (A.14) 
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r 

^ggss 



^ - ' (5 - k)x^x^o~ + (5 + /i)x+X"c++ (A.15) 



^-2ibsw'X^hc^ + 2ksvfX^X3C'^ — ^/i/ic"*" — ft^XsXsc^ j + h.c 
+ f-^X^^c" + «X"^c+ + x"^X3C" + X~X3C+ + — XaXsC^ - —hhc^ 

B Interaction vertices 

We give the full list of the propagators and interaction vertices for the Standard Model 
with non-linear gauge-fixing. These results follow straightforwardly from the Lagrangian 
written above: the coefficient in front of the fields is multiplied by i and we do a transfor- 
mation to momentum space by substituting 9^ Also there comes a factor from 
permuting identical fields. 

This list can be checked with the one given in the Appendix of [5J. The non-linear 
gauge-fixing used there is the same, and so the resulting terms are the same. The only 
difference in conventions there is that the Goldstone bosons x^^Xs ^ire defined with an 
additional minus sign, so all the interaction vertices that involve an odd number of Gold- 
stones have a relative minus sign. Also note that |[5j has the momentum of the antighosts 
defined with an opposite sign and does not include the overall i factor for vertices. 

Notes about fermions: some vertices are given generally for any fermion / in terms of 
its charges. The charges used are the charge of the left-handed component: 

t%UL,i^e) = +l, t=^K,ei) = -^, (B.l) 

and the electric charge q: 

q{u) = ^, qid) = -^, qii^e) = 0, g(e) = -1. (B.2) 

All fermion vertices apply to all 3 generations: 

e^(e,/i, r), z/g ^ (z/e, i/^,, z/^), u^{u,c,t), d^{d,s,b), (B.3) 

with the same charges, but different masses. As it was mentioned, SU (3) is not discussed 
in this work, but note that for each quark there are 3 states of different color. Also note 
that (1 — 7^)/2 and (1 + 7^)/2 are the projectors of the left-handed and right-handed 
components, respectively. 
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Propagators 

(X^X^) = ^^3^' (B.5) 
= (B.6) 

(//) = (B.8) 

'g — rrif 

Vector- Vector- Vector- Vector 

W^W-A,A, ie' {-29^" g''' + (1 - a'/^w) {9"' 9"^ + 9'" 9"')) 

W^W-A,Z, ie^^ (-2^'^'^/'^ + (1 - a~^liw){9^'9'"' + ^'^'^/O) 

2 

ie^ Y {-29'''' 9'" + (1 - ^Vei^)(^'^V + 



Vector- Vector- Vector 

'^;{Vx)Wt{ji2)A,{vz) le \9^''{Vx - ^2)^ + 9'''{V2-Vzr + ^"^(Pa -Pi)'^ 

W^;(^^i)^^^(^^2)^p(p3) \9^\vx - + 9'''{P2 - PsT + 9"'{P3 - PiT 

Sw 



+ ~PKwW-P29 



Vector- Vector-Scaleir-Scalar 
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Z^iA^X'^X^ 

z^z^x'^x' 



w+w-hh 



W^W-X3X3 

w+w-x+x- 



Z^Z^hh 



Z^Z^XsXs 



-ie^- — (1 - ak)g 



jJLV 



1 



±e' — {l-f55cl,/s^^)g^^ 
2ie''g^"' 

CwSw 



2c]ySw 



1 



1 



1 



1 



1 



le 



Vector- Vector-Scalair 



ZfxZyh 



^emw{^ — 0()g' 



Cw 

le g 



le 



Sw 

2 9 



Ve ct or- S caleir- S calEir 



X~(Pl)X^(P2)A^(p3) 

X~{Pl)X^{P2)Zf,{p3) 

h(Pl)X3(P2)Z„{p3) 



e^{p'2-p1 + h's) 
±ie-^{p^ - p'l + Rp'^) 
«e(p2 -Pi) 



S calar- S calar- S calair- S calair 



hhhh 

X3X3X3X3 

X±X±X=fX=f 
hhxsXs 
hhx^X~ 



X'^X X3X3 



—le 



-le 



2 J3m|_ 



m 



-le 



H 



—le 



2 

2 

rrii 



1 + 2e'izm'z/m'H) 



-^e" . , (1 + 25^^vK?7i^/m^) 



m 



H 



1 + 2k Cw^m^/m^) 



S caleir- S calar- S caleir 



hhh 

hx'X'^ 

hX3X3 



-le- 



m 



-le- 



H 



2sw'>Tiw 



m 



—le- 



H 



[I + 25iwrrv^w/m]j) 



Fermion-Fermion- Vector 



IfA, 
ffZ, 

i7eeW+, eueW-, udW+, duW 



le 



SwCw 

1 

\/2sw 



^{l - 7^)/2 



Fermion-Fermion-Scalcir 
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ffh 

eexs, ddxs 

UUX3 

udx" 
dux~ 



nif 



— 26 " 

5 

2swfnw^ 



+e- 



\/2sw'>Tiw 



(l + 7')/2 
■(l + 7^)/2 



y/2sw'mw 

-e^^ {ma{l + 7') - - 7^)) /2 



\/2swmw 
1 



(m„(l + 75)-md(l-75))/2 



Ghost-Ghost- Vector 



c'^{piy{P2)w^ 

c'^{pi)c^{p2)A^ 



Sw 

±ie{p^ + ap^) 
±ie'^{pi + M) 

Sw 

Tie{p1 - ap^) 

Sw 



Ghost- Ghost-S calEir 



c^c^h 



2swCw 
Czmz 

2sw 
^emwiw 
mwi 



—le 



±e 



(1 + e) 



w , 2 



±e 



2swCw 
. mwiw 

2sw 
mwiw 



-le- 



(1 + ^) 



2sw 



{l-h) 



Ghost- Ghost- Vector- Vector 



Sw 



c'^c^A^A^ 2ie^ag^"' 



,2 



Sw 



c-^c^Z^Z, 2ie^~P^g^'' 



Sw 
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Ghost- Ghost- S calEir- S calai 



c^c^hh -ie^e^ 



-Z Z . 2~ 

c c X3X3 «e £ 



2s2 r2 



c c X+X3 «e e—^ 

-=F A ± -2- ^H/ 

c+c X X3 «e 5 + K{cy^ - s^) 1 



-=F ±1,1, -2? 

z^c hh —le 0—^ 



-=n ± • 2~ 'Cm/ 



c^c X3X3 -«e ^2^2 



vv 



-=F±i, _l2/J ~\ 'Ch^ 

c c X3'i ±e (o — 



4s 



2 



c+c X X «e (d + Ac)^ 



-=F =F ± ± • 2/r ~\ 

c+c+x X «e 



C FeynArts model files 
C.l Lorentz.gen 

Here we list the modified pieces of tiie Lorentz . gen generic model file. These cor- 
respond to Vector- Vector- Vector, Vector-Scalar-Scalar, Ghost-Ghost- Vector- Vector and 
Ghost-Ghost-Scalar-Scalar vertices respectively. 

(* V-V-V: *) 



AnalyticalCouplingE si V[jl, moml, {lil}] , s2 V[j2, mom2, {li2}] , 
s3 V[j3, mom3, {liS}] ] == 
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G[-l][sl V[jl], s2 V[j2], s3 V[j3]] . 

{ MetricTensor [lil , li2] FourVector [mom2 - moml, li3] + 
MetricTensor [li2 , liS] FourVector [momS - mom2, lil] + 
MetricTensor [113 , 111] FourVector [moml - mom3, 112], 
FourVector [moml , lil] MetricTensor [112, 113], 
FourVector [mom2, 112] MetricTensor [113, lil], 
FourVector [mom3, 113] MetricTensor [lil , 112] } 

(* S-S-V: *) 

AnalytlcalCoupllng [ si S[jl, moml], s2 S[j2, mom2] , 
s3 V[j3, momS, -[113}] ] == 
G[-l][sl S[jl], s2 S[j2], s3 V[j3]] . 
{ FourVector [moml - mom2, 113] , 
FourVector [mom3, 113] } 

(* U-U-V-V: *) 

AnalytlcalCoupllng [ si U[jl, moml], s2 U[j2, mom2] , 
s3 V[j3, momS, {113}], s4 V[j4, mom4, {114}] ] == 
G[l][sl U[jl], s2 U[j2], s3 V[j3], s4 V[j4]] . 
{ MetricTensor [113, 114] } 

(* U-U-S-S: *) 

AnalytlcalCoupllng [ si U[jl, moml], s2 U[j2, mom2] , 
s3 S[j3, mom3] , s4 S[j4, mom4] ] == 
G[l][sl U[jl], s2 U[j2], s3 S[j3], s4 S[j4]] . 
{ 1 } 



C.2 SlVE.mod 

We give some examples of the modified vertices as defined in the classes model file. 
Note the new terms proportional to the gauge-fixing parameters that are called Galpha, 
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Gbeta, Gdelta, Gkappa, Gepsilon. The (W+W-AA), {W+WA) and {W+Ax~) ver- 
tices that were needed in the tree-level 77 — > W'^W~ calculation are defined as: 

C[ -V[3], V[3], V[l], V[l] ] == -I EL^2 * 

{ {2, 4 dZel + 2 dZWl + 2 dZAAl - 2 CW/SW dZZAl}, 

{-1 + Galpha^2/GaugeXi [W] , -2 dZel - dZWl - dZAAl + CW/SW dZZAl}, 

{-1 + Galpha~2/GaugeXi[W] , -2 dZel - dZWl - dZAAl + CW/SW dZZAl} } 

C[ V[l] , -V[3] , V[3] ] == -I EL * 

{ {1, dZel + dZWl + dZAAl/2 - CW/SW dZZAl/2}, 
{0, 0}, 

{- Galpha/GaugeXi [W] , 0}, 
{+ Galpha/GaugeXi [W] , 0} } 

C[ S[3], -V[3], V[l] ] == -I EL MW * 

{ {1 - Galpha, dZel + dMWsql/(2 MW"2) + dZWl/2 + dZAAl/2 + dZGpl/2 + 
SW/CW dZZAl/2} } 

The coefficients of the interaction vertices are listed in the first column, while the second 
column involving renormalization constants dZel, dZWl, etc. lists the coefficients of the 
counterterm vertices, which we have not modified. Following are examples of other vertices 
that had new kinematic terms introduced, which arc Vcctor-Scalar-Scalar and 4-vertices 
with ghosts. We give here {hx~W~^), {c~c^ZW~^) and (c^c^hh): 

C[ S[3], S[l], -V[3] ] == -I EL/(2 SW) * 

{ {1, dZel - dSWl/SW + dZWl/2 + dZHl/2 + dZGpl/2}, 
{Gdelta, 0} } 

C[ -U[4], U[l], V[2], -V[3] ] == I EL~2 (CW/SW) /Sqrt [GaugeXi [W]] * 
{ {Gbeta} } 

C[ -U[2], U[2], S[l], S[l] ] == -I EL~2 Sqrt [GaugeXi [Z] ]/ (2 SW^2 CW~2) * 
{ {Gepsilon} } 
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